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, Abstract 

, Supersymmetric method of the constructing well-like quasi exactly 

solvable (QES) potentials with three known eigenstates has been ex- 
tended to the case of periodic potentials. The explicit examples are 
. presented. New QES potential with two known eigenstates has been 

' obtained. 



1 Introduction 

Description of the electron's motion on a lattice has been investigated for a 
long time as a central problem of the condensed matter physics. Such quan- 
tum problem is reduced to the solving of the Schrodinger equation with some 
model potential which is periodic often. Therefore, the periodic quantum me- 
chanics problems remain at the investigation's focus up to now. 

The general properties of the solutions of Schrodinger equation with pe- 
riodic potential energy are described by the oscillation theorem ^J. Energy 
spectrum of the periodic potential has band structure, i.e. eigenvalues belong 
to the allowed bands (energy bands) [Eo,E{\,[E^,E2],.... The wave functions 
are the Bloch functions, which are bounded and extended on the full real axe 

ipix + L) =exp'''^^{x), (1) 
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where L is potential period and is a so-called quasi-momentum. The limits 
of the energy bands are given by the equation kL = {0,7?}, and the wave 
functions, which belong to the limiting energy values, satisfy the condition 
ip{x + L) = ±tlj{x). These energy values and wave functions are often called 
the eigenvalues and the eigenfunctions of the described above problem. 

Oscillation theorem claims that in the case of periodic potentials the 
eigenfunctions, which belong to the limits of the energy bands and are ar- 
ranged in the energy of the increasing order Eq < Ei < Ey < E2 < E21 < 
E^..., are the periodic functions with the period L, 2L, 2L, L, L, 2L, 2L, ... and 
have 0, 1, 1, 2, 2, 3, 3, ... nodes in the interval L respectively. 

Despite long term investigations, there is rather a limited number of ex- 
actly solvable periodic potentials even in one dimension. The classical exam- 
ples are the Kronig- Penney model potential *?] or Lame's potentials pi. 

Because of limited number of the exactly solvable potentials, recently 
much attention has been given to the quasi exactly solvable (QES) poten- 
tials for which a finite number of the energy levels and the corresponding 
wave functions are known explicitly. A general treatment of the quasi exact 
solvability has been introduced by Turbiner and Ushveridze The class of 
QES trigonometric potentials was presented in In jH] it was shown that 
the Lame equation is a peculiar example of QES systems. The authors of 
the paper [7j considered a family of spectral equation which extends those 
of [U]. Authors of jH] has applied quantum Hamilton- Jacobi formalism to the 
QES periodic potentials. In the latest paper |S] an unified treatment of quasi 
exactly solvable potentials was proposed. 

The powerful tool for studying the problem of exact solvability of the 
Schrodinger equation is the supersymmetric (SUSY) quantum mechanic in- 
troduced by Witten (TU] (for a review of SUSY quantum mechanics see [TT]). 
The SUSY method for constructing QES potentials was used for the first 
time in jT^j - fl] . The idea of this method starts from some initial QES po- 
tential with n + 1 known eigenstates and using the properties of the unbroken 
supersymmetry to obtain the SUSY partner potential, which is a new QES 
one with n know eigenstates. 

In |15j-|r7j using the formalism of SUSY quantum mechanics a large num- 
ber of new solvable and QES periodic potential was proposed. It is worth 
mentioning recent paper JH] , where the highest order SUSY transformations 
was applied for studying periodic potential. 

In recent Tkachuk's papers 1201-1221 a new SUSY method for constructing 
of the QES potentials with two and three known eigenstates has been pro- 
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posed. This method does not require knowledge of the initial QES potential 
in order to generate a new QES one. Within the frame of this method QES 
potentials has been obtained for which the explicit form of the energy levels 
and the wave functions of the ground and the excited states can be found. 
After the paper [221 Dolya and Zaslavskii, where they showed how to gen- 
erate QES potentials with arbitrary two known eigenstates without resorting 
to the SUSY quantum mechanics, SUSY method has been extended j21| for 
constructing QES potentials with arbitrary two known eigenstates. In our 
recent works using the SUSY method periodic [25j and disordered [2^] QES 
potentials were obtained. 

In the present paper using the results of previous study [201- EH] we extend 
Tkachuk's SUSY method for constructing QES periodic potentials with three 
known eigenstates. 

2 The Witten model of SUSY quantum me- 
chanics 

Witten model of supersymmetric quantum mechanics is a quantum mechan- 
ics of the matrix Hamiltonian 

(2) 

where Hamiltonians 

H^ = -l^^+V^{x)=B^B^ (3) 
are supersymmetric partners and 

Here h = m = 1 units are used. Function W{x) is referred to as superpoten- 
tial, V±{x) are so-called supersymmetric partner potentials 

2V±{x) = W\x)±W'{x). (5) 

Energy spectrum of the supersymmetric partners and is identical 
except for zero-energy ground state which exists in the case of the unbroken 
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supersymmetry. This leads to twofold degeneracy of the energy spectrum of 
if, except for the unique zero-energy ground state. Only one of the Hamil- 
tonians H± has zero-energy eigenvalue. We shall use the convention that the 
zero-energy eigenstate belongs to if_ 

=0 ' 

where n = 0, 1, 2, .... The wave functions of the supersjTiimetric partners H± 
are related by the supersymmetric transformations 



V ^n + l 



(7) 



Due to the factorization H_ = B^B , we can find solution of the Schrodinger 
problem for the eigenstate with zero energy 



It is easy to see that 



H.tlj^{x) = Eo^oix)=0- (8) 



ipQ (x) = Cq exp (- J W{x)dx), (9) 



where Cq is an arbitrary constant. 

In the present paper we shall consider the systems on the full real axe 
— oo < X < oo with periodic superpotential. The periodic superpotential 
W{x + L) = W{x) leads to the periodic potential energy V±{x + L) = V±, 
which results in the bounded and extended eigenfunction. A satisfactory 
condition for the existence of periodic eigenf unctions, written in the terms of 
the SUSY quantum mechanics, is 

L 

W{x) = 0. (10) 







In 123 12H1 a detailed analysis of the SUSY quantum mechanics was made for 
this case. 
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3 SUSY constructing QES potentials 



We shall study the Hamiltonian iJ_ with the potential energy 

V,{x)=W^,{x)-W',{x), (11) 

the ground state of which is given by ©• 

Let us consider Hamiltonian which is the SUSY partner of Hamilto- 
nian . If we calculate the ground state of if+ we immediately find the first 
excited state of using the degeneracy of the spectrum of SUSY Hamilto- 
nian and SUSY transformations ((Tj). In order to calculate the ground state 
of let us rewrite Hamiltonian in the following form 



where 



= + e, e > (12) 
= BiB^, (13) 

and Wi{x) is a some new function. Note that e is the energy of the ground 
state of since has zero-energy ground state. 

The ground state wave function of if+ with the energy E = e is also zero 
energy wavefunction of and it satisfies the equation 

B^^+{x) = 0. (15) 

The solution of this equation is 



V^+(x) 



C+ exp - y Wi (x) dx^ , (16) 



where Cq is an arbitrary constant. 

Using the SUSY transformation (|7j) we can calculate the wavefunction of 
the first excited state of H_ . Repeating the described procedure for we 
can obtain the second excited state for H_ and so on. This procedure is well 
known in the SUSY quantum mechanics (see review for example). The 
wavefunctions and corresponding energy levels read 



i^ni^) = C^B+...B+_2B+_^exp(^- fWn{x)dx 



(17) 
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where n = 1,2, N; eo = e, Bq = B^, Wq{x) = W{x), C~ are an arbitrary 
constants. Operators B^ are given by (jD) with the superpotentials Wn{x). 
Equation ()12p rewritten for N steps 

i/W + e„, (18) 

where n = 0,1, N — 1, leads to the set of equations for superpotentials 

W^{x) + C(a^) = W^n+i(^) - + 2en, (19) 

where ?t, = 0, 1,...,A^ — 1. 

Unfortunately, each of the equations in (fT^ are the Rikatti equation, 
which can not be solved in the general case. Previously this set of equations 
was solved in special cases of shape-invariant potentials [21] and self-similar 
potentials for arbitrary (see review jHO])- For = 1 in the context of 
parasupersymmetric quantum mechanics one can obtain a general solution 
of ()19|) without restricting ourselves to shape-invariant and self-similar po- 
tentials 31J. In recent papers 1201-1211 a solution of (fT^ for = 1 and N = 2 
in order to obtain non-singular QES potentials with two and three known 
eigenstates respectively has been constructed. 

Let us write set of equations (jl9j) for the case = 2 in the explicit form 



W^{x) + W^(x) = W^{x) - W[{x) + 2eo 
Wf{x) + W[{x) = wf{x) -Wf,{x) +2ei ' 

It is convenient to introduce new functions 

W+{x) = Wi{x) + Wo{x) ( W+{x) = W2{x) + Wi{x) 
W_ lx) = Wi {x) - Wolx) \ W^ {x) = W2{x) - Wi {x) 



(20) 



(21) 



then superpotentials can be rewritten in the following form 

(22) 

In the terms of new functions ()21|) the set of equations (j^n|) read as follows 

Wiix)=W^{x)W+ix) + 2eo . . 

W\_{x)=W_{x)W+{x) + 2ei ■ ^ ^ 



2Wo{x) = W+{x) - W^{x) ( 2Wi{x) = W+{x) -W.{x) 
2Wi{x) = W+lx) + W^lx) \ 2W2{x) = W+{x) + W^x) 



6 



Note, that there are two terms for the Wi{x) in the equations (|22j) with 
respect to W^{x) and with respect to W+{x). This gives us a possibihty to 
obtain relation between W+{x) and W+{x) 

, , W'{x)-2e^ ~ , , W'{x)-2e^ , , 

«'-w^^k^-^''-w-T4r' '''' 

here (|^ are used. It is easy to rewrite this equation as follows 

W^{x)W^{x)^^{x) - W^{x)\ - [l^+(a:)iy+(x)]'+ 

+2[eiiy+(x) + eoW^+(a;)] = 0, 



or 

where we have introduced a new function 

U{x) = W+{x)W+{x). (26) 

We arrive again to the Riccati equation with respect to U{x). On the other 
hand, this is an algebraic equation with respect to W+{x), which can be 
solved explicitly 

W - 2l/(x)(^(x)+2eo) 

+ U'(x){l+R(x)) /^yN 

/ \ _ U'{x){l+R{x)) > 



2(U(x)+2eo) 



where 



f/(x)(f/(x)+2eo)(t/(x)-2ei) 



lH(x) = l + 4 "^^^^"^^^ ^' ^ ^(x) = ±V^- (28) 

The square root yi{x) is a positively defined value, while the function R{x) 
can be chosen in the form of lH(x) or — 9^(x) within different intervals sepa- 
rated by zeros of the function d\{x). 

Thus, we can start from an arbitrary function U{x) to construct the func- 
tions W+{x) and W^{x) given by Using (j221) we obtain three consequent 
superpotentials 

Woix) = Uw4-)-^^^ 

W,{x) = l{w,{x) + ^^) . (29) 



W+{x) 



Then because of (jl7p . we can find the wavefunctions of three exphcitly known 
eigenstates of the Hamiltonian H_ 



^r(x) = Cfiy+(a;)e-/'^i(^)'^^ (30) 
-ilj-(x) = ({Wo{x) + W2{x))W+{x) - iy|(x))e-/^2Wrf^ 

where energy values are Eq = 0, = eo, E^ = eo + ei and potential energy 

V4x) = \{Wo{xr-W;,{x)). (31) 

Simultaneously we can find the wave function of two explicitly known 
eigenstates of the Hamiltonian 

^tix) = B^ij^ix) ^^^^ 
with energy values E^ = eo, E2 = eo + ei and potential energy 

V.ix) = ^iWoixy + W;^ix)). (33) 

Note that obtained terms for the superpotentials, potentials and wave 
functions allow existence of two different solutions depending on the selected 
sign before the square root ±^/R{x) in the W+{x) and W+{x) definitions. 
Here and later we shall distinguish solutions which were obtained for different 
signs, by superscript in the parenthesis after the function designation, for 
example Y(x)^^\ We will denote as Y(x) solutions which are identical for 
different signs Y{x)^^^ = Y{x)^^\ 

Choosing different generating functions U{x) we will obtain different QES 
potentials (|31|) with three explicitly known eigenstates (|3Up and QES poten- 
tials with two explicitly known eigenstates ((221) • Of course, function U{x) 
must satisfy some conditions to provide physical solutions of the Schrodinger 
equation. 

The main obvious condition imposed on the function U{x) is a positivity 
of the expression under the square root 

. , 4[/(a:)(^(a;)+2eo)(^(a;)-2eO 

U'(xY - 
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on the all periodicity interval. 

Another set of restrictions imposed on function U (x) appears due to the 
requirement of the non- singularity of resulting potential V-{x). The full anal- 
ysis of the properties of superpotential Wo{x) which provides non-singular 
potential V-{x) and wave functions ipQ^x), ipi{x) in the terms of function 
W+{x) was done in [211- (SHI for the case of the quasi exactly solvable po- 
tentials with two exactly known eigenstates. Below we extend this analysis 
to the case of the quasi exactly solvable potential with three exactly known 
eigenstates. 

As we can see from the superpotentials Wo{x), Wi{x) and W2{x) defini- 
tions, potential V_{x) can have poles at the points Xq where W+{xo) = or 
W^{xo) = 0. Fortunately, such poles can be removed when [21] 

W\_ixo) = ±2eo, , . 

iy;(xo) = ±2ei. ^""^^ 

Besides, potential energy V- (x) can have poles at the points of singularity 
Xoo of the function W^{x). As it was shown in pOlj, if function W^{x) at the 
singularity points x^o has the behavior 

W+{x) = const H h o{x — Xoo), (36) 

X .^-oo 

or 

_3 

W+{x) = — + o{x-x^), (37) 

X Xqq 

obtained potential energy and wave functions will be continuous functions at 
the points x^o- 

To provide bounded and extended wave functions ipQ^x), ipi{x), ^/'^(a;) 
the conditions (fTUj) should be satisfied 

j^WQ{x)dx = 

J^Wi{x)dx = . (38) 
J^W2{x)dx = 

These conditions are satisfied in the simplest way if the corresponding super- 
potenials Wo{x), Wi{x), 1^2(2;) are odd function with regard to the middle 
of the periodicity interval Xm- To obtain odd superpotentials it is enough to 
expect the odd behavior of the function W+{x). 
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Let us choose the U{x) as even function with regard to the middle of the 
periodicity interval Xm- Then, if we apply solutions with the different signs 
before square root to the parts of the periodicity interval from the left and 
from the right of Xm, W+{x) will be odd function. It is easy to see if we 
rewrite term ^I7\ for the W+{x) as follows 

U'{x) ± ^/U'{xy + 4U{x){U{x) + 2eo)(f/(x) - 2ei) 

Application of the solutions with the different signs leads to the finite 
breaks of the function W+{x). These breaks can be removed if the value of 
the function W^{x) will tend to zero both from the left and right direction. 

Thus, to provide existence of the bounded extended wave functions, U{x) 
should be even function with regard to the middle of the periodicity interval 
Xm, and obtained function W+{x) should have zero at the point Xm- Function 
can have zeros at the points, where U{x) = 0, and, since U{x) should 
be even function with regard to Xm, function U{x) can have at the point Xm 
zero of the even-order only. 

Of course, function U{x) can have zeros at the other points of the periodic- 
ity interval too. Let us analyze in details the behavior of the superpotentials, 
potential energy and the wave functions in the vicinity of the U{x) zeros. 

Let the function U{x) have the first-order zeros at the points Xq 

U{x) = U'{xl){x - xl) + \u"{xl){x - xlf + o{x - (40) 

Then behavior of the functions Wj^{x\ Wj^{x) in the vicinity of the points 
will be as follows 



(41) 



It is easy to see that functions iy+(x) and W^^ix) at the points Xq will have 
non-zero values or will have zeros which satisfy 

Superpotentials Wq{x^, W\{x), W^ix) will be the following 



(42) 







(±) 


- ^0 


+ o{x — 


Xq 




K-^ , 


i±) 


= ) 


+ o{x — 


Xq 


W2 


[< 


(±) 


= 4") 


+ o{x — 


Xg 
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where 

Ai+) _ Ai+) _ 8egei+C/'(a:g)2-eoC/"(xg) 

^1 - ^2 - 2eii7'(i§) ■ y^"^) 

Ai-) - .(+) _ U"(x^)-8eoei 

^0 — ^2 — 2!7'(xg) 

Obtained potential will be regular function too 

V_ {x) = ai^^ + o(a; - x;;) , (44) 



where 



a 



(-) _ 64c^q+8ri;;\x;;)-+;;''(x;;)--i(ko((;'(x;;)-+<:ir''(x;;))-2f;'(,^^^^ 



= -3ai-) + 460 + 



8;7'(a;g)2 



2C/'(a:g) 



(45) 



The wave functions ipQ (x), ipi (x), ip2 {x) will read as follows 

i>Q{x) = 1 + 0{x - Xq) 

^-{xY+) = 2eo{x - x^) + o{x - x^y 
^-(a;)(-) = £^ + o(x-a;g) ■ ^^^^ 

ip2{x) = — 2ei + o{x — Xq) 

Thus, at the points x^, where function U{x) has first-order zeros, potential 
V^{x) and wave functions iIjq{x), ipilx), ip2{x) will be continuous functions, 
and wave function '4'i{x) can have simple zeros at the points Xq depending 
on the selected sign before the square root. 

Let us consider potential V+{x), which is the supersymmetric partner of 
the obtained potential V-{x). Potential V+{x) will be regular function in the 
vicinity of the points Xq too 



T/+(x)W = Q;f^ + o(x-xS) 

„(+) _ ^,(-) ^ ^ C/'(xg)2-2.oC/"(a:g) 



ay ' = ai-) + 261 + ^'(-o)--2.o^-(a.g) (47) 



a\ ^ = q;L+^ - 2ei 



with continuous wave functions 



'tpi{x) = A/2eo + o{x - Xq) 

^2+(a:)W = 2v^ei(eo + ei)(a;-a;g) + o(a;-a;g)2 . (48) 
^+(^)(-) = (-°+^^J("g) + o{x - xg) 
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Depending on the selected sign before the square root wave function 
can have nodes at the points Xq. 

Now let the function U{x) have second-order zeros at the points Xq 

U{x) = lu"{xl){x - xiy + lu^'\x'^){x - xy + o{x - xl)\ (49) 

then behavior of the functions W+{x) and W+{x) will be as follows 

W+ix) = 26o(x - 4) + oix - 4)3/2 
W+ (x) = 2ei {x-xl) + o{x - x[5)3/2 ' ^"""^ 

i.e. at the points Xq functions W^{x) and W+{x) will have zeros. Keeping in 
mind ()35|). it is easy to obtain the following coefficient restriction 

= Wl{xl)W+{xl) + 2Wi{x'f,)W[{xl) + W+{x'o)W:i{xl) (51) 
= 2M/;(4)iy;(x^) =8eoei. 

Note, that existence of the fractional powers in the series expansion leads 
to undesired poles of Wo{x) at the points Xq 



^''»<^''^'-4v'SI)+"'^-^»'"- '''' 

which can bring the singularity to the potential energy V- (x). It is easy to see, 
that in the case of U^^\x^o) = the fractional powers in the series expansions 
disappear 

W+{x) = 2eo{x - x'^o) + o{x - x'^of 
W+{x) = 2ei(x-x[;) +o(x-x[;)2 ' ^ ^ 

The condition U^^\xq) = is satisfied in the simplest way if the point Xq is 
a middle of the periodicity interval and U{x) is even function with regard to 
the Xq, which at the same time provides the fulfilment of 
Then superpotentials read 





[< 


)'^) 


- B^^^ 

— -°o 


+ o{x — 








)'^) 


= 


+ o{x — 


Xq 


W2\ 




)'^) 


= 5f ) 


+ o{x — 





(54) 
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where 

d(+) _ d(-) _ r(+) _ D 

B^-'^ = = B^-'^ = -B , (55) 

B = l/4v/32(6o-6i) + f/(4)(,T[5)/(2eoei) 

Obtained potential V-{x) will be continuous function 

= /?W + o(x-a;S) 

64eoei + 320€o€iS 

Wave functions ipo{x), '^i{x), 1/^2 {^) will read as follows 

i/joix) = 1 + o{x-Xq) 

V'r(x) = 2eo(x - x^) + o(x - x^)2 . (57) 
V'^(x) = -2ei + o{x - x^q) 

Thus, in the vicinity of the second-order zero of the U (x) potential energy 
V-{x) and the wave functions ip^lx), ipi{x), ip2i^) ^^^^ be continuous func- 
tions, if Xq — Xjn is the middle of the periodicity interval and U (x) is even 
function with respect to Xq. Wave function ipi{x) will have node at the points 

The supersymmetric partner V^{x) of the V-{x) potential in the vicinity 
of the Xq will have the following behavior 

V^{x)(^^ = P^^^ + o{x - x'q) 

P+ ~ ^0 ^tl 64eoei ^ 320eoeiB 



with the following wave functions 



i)t{x) = V2eo + o{x-xl) 

ip+{x) = 2V2ei{eo + ei){x-xl) + o{x-x^Q] 



2 



(59) 



Thus, in the vicinity of the second-order zeros Xq of the function U (x) poten- 
tial energy V+{x) and the corresponding wave functions iljf{x), ip2{^) will be 
continuous function and wave function '0^(x) will have nodes at the points 



Xn- 



Let us analyze the case when the function U{x) has the highest order of 
zeros at the points Xq using the particular case of the third-order zeros 
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U{x) = lu^'\x:){x - xlf + ^^U^'\x:){x - xlf + o{x - xlf. (60) 

Then the series expansion for the function W+ {x) in the vicinity of the points 
Xq will start from the terms which will be proportional to the (x — Xq)'^''^, 
thus, condition ()35|) will not be satisfied, and then obtained potential energy 
V_{x) will have poles at the points Xq. Consequently, function f/(x) should 
not have zeros of the highest then second orders. 

Singularities at the potential energy, except the zeros of U{x), can appear 
at the points where U'{x) = or 1 — \jR{x) = 0, that is 



U'{x) 


= 0, 


U{x) 


= 0, 


U{x) 


= -2e 


U{x) 


= 2ei. 



0, 



(61) 



Case of U{x) = was considered in the details above. In the vicinity of 
the points Oq, where the derivative of f/(x) is equal to zero, i.e. U'{ao) = 
and f/(ao) 7^ 0, generating function U{x) can be written as 

U{x) = U{ao) + if/"(ao)(x - aof + o(x - aof. (62) 
Then behavior of function W+{x) in the vicinity of ao will be the following 



. ^(±) , / ^(ao)(2eo + [/(ao)) Wjao) ^ ^ 

^■^("^ = Uiao) - 2e, + 4e,-2Uia/ ' ' + " ' 

(63) 

in other words, in the vicinity of zeros of U'{x), which do not coincide with 
zeros of U{x), obtained solutions will be continuous functions. 

In the vicinity of bo, where U{bo) = 2ei, function W+{x) will behave as 
follows 

w^i-y^^ = ^^+o{x-bo), 

iy+(x)(-) = + const + o{x- bo). 

Despite singularity of the function W+{x)^^\ potential energy and wave func- 
tions will be continuous functions, because pole of W+{x) satisfies the con- 
dition (jnni)- 
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In the vicinity of Cq, where U{co) — — 2eo, function W+{x) can be figured 
out as foUows 



thus, potential energy and wave functions will be continuous functions again. 
Consequently, at the all points, where denominator of W+{x) can turn into 
zero, potential energy V-{x) and wave functions iPq{x), ipi{x), ■02" (^) ^^^^ be 
free of singularities. 

Similar analysis, which we shall omit due to its inconvenience, with re- 
spect to the potential V+{x) shows, that potential V+{x) and corresponding 
wave functions will be free of singularities at the all considered points ex- 
cept the points cq, where potential V+{xY'^^ will have pole with the following 
behavior ^ 

V-^.ix)^'^^ — -. r- -|- const -\- oix — Co). (66) 

[X - CqY 

Fortunately, this singularity can be avoided if within the parts of periodicity 
interval which contains co we apply solution V^{x)^~^ instead of Vj^{x)^'^\ 

Another way to avoid singularities in the potential V+{x) is to exclude ze- 
ros in the denominator of W+{x) by picking up the amplitude of the function 
U {x) in such a manner that equations U {x) + 2eo = and U {x) — 2ei = 
not be fulfilled. Indeed, since energy levels eo, Ci are positively defined values 
and U {x) is a periodic bounded function, we can always fit the amplitude of 
generating function U {x) using the following rule 

eo < l/2min[/(a;), , . 

ei > l/2maxC/(x), ^ ' 

where min U {x) and max U {x) - minimal and maximal values of the U (x) at 
the periodicity interval respectively. 

Thus, periodic function U{x) generates quasi exactly solvable potential 
V-{x) with three known cigenfunctions i/jqIx), il)i{x), ip2{x) for the energy 
values Cq > and ei > 0, if R{x) > for all periodicity interval. Simulta- 
neously, function U{x) generates quasi exactly solvable potential V-^{x) with 
two known eigenfunctions ipti^), ^^(a;) in the case of U{x) G (— 2eo; 2ei) 
and R{x) > for all periodicity interval. 

To provide free of singularities potential energy and extended bounded 
wave functions, U{x) must be even function with respect to the middle of 
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the periodicity interval Xm and must have second order zero at this point. 
Generating function U{x) may have first-order zeros at the other points of 
the periodicity interval and should not have zeros of the highest order. The 
derivative of the U"{x) at the point Xm should satisfy the condition U"{xm) = 
SeoCi. It is necessary to use solutions with opposite signs from the left and 
right sides with regard to point Xm- 

To illustrate the above described method we give a short example. 

Trigonometric extension of the Razavy potential. Let us start from 
the generating function 

U{x) = 4eoei sin^ x. (68) 

Similar generating function U{x) = 4eoei sinh^ x at ei = eo + 1/2 gives well 
known quasi exactly solvable Razavy potential [22] ■ Than can be rewrit- 
ten in the following form 

«H(x) = (-1 + 2eo - 2ei + 4eoei sin^ x) tan^ x. (69) 

We shall omit the general expression for the superpotentials and potential 
energy as it is huge and rather useless. There are at least three sets of eo, 
ei, which allow us to resolve the root in the function R{x) and therefore to 
significantly simplify the final results. 

The first set is 

4eoei = 
-1 + 2eo - 2ei > ' 



(70) 



for which we obtain trivial solution eo = or ci = 0, what leads to the 
U{x) = 0. 

In the case of the second set 



-l + 2eo-2ei = -4eoei 
-1 + 2eo - 2ei > 



(71) 



we obtain ei = —1/2. Then 



WA.) = ''p'. . (72) 

1 + V 2eo sin X 

Function W^{x) has zeros at the points Xk = 7m/2, n = 0, ±1, .... The deriva- 
tions W!^{x) at this points are — 2eo/(l + i/2eo) or 2eo and condition pSj) is 
not fulfilled. 
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The last set 



-l+2eo-2ei = 
4eoei > 



(73) 



gives ei = eg — 1/2, then square root can be rewritten in the following form 
Function W+(x) reads as follows 



R{x) = 2^eoei sin x tan x 
eo > 1/2 



W+{x) 



2eo(cos^ X + 2eo sin^ x) tan x 



1 + 2y/eoei sin x tan x 



(75) 



Function W+(x) has singularities at the points x[.^-* = iarccos ^/eo/ei + 
2Tm, n = 0, ±1, ... and x^^^ = ± arccos(— A/ei/eg) +27rn, n = 0, ±1, .... Due to 
the limitation eo > 1/2, solutions x^,^"* belong to the complex space and thus, 

(2) 

can be dismissed. At the points x^ function iy+(x) has simple poles with 
the pole coefficient —1, thus potential energy VL(x) will be regular function 
at points x^^^ for any eo. Additionally, function W^{x) has simple zeros at the 
points Xfc = vrn, n = 0, ±1, .... The derivations W+{x) at all these points are 
equal to 2eo, so all conditions imposed on generating function U (x) to provide 
non-singular real potential energy V^{x) are satisfied for any eo > 1/2. 

Then, using the definition of function W+{x) (jTTj) . solution for super- 
potentials Wo{x), Wi{x), W2{x) and relation between superpotential 
Wo{x) and potential energy V-{x) (fTTjl . we can find three eigenstates of the 
potential 



(x) = eo - ^ + ^ (eoei - 6 v^e^ cos x - eoei cos 2x^ , 

where ei = eo — 1/2. The energy values of this eigenstates are Eq 
= eo, E2 = eo + ei and wave functions are given by 



(76) 



ijoix) = Co-e^^=°^'i(l + 4(yi^ + ei)cos2f) 
^-(x) = Cfe^'^^^'feosinx . (77) 

^^{x) = C2-e^^™«'i2ei(l+4(yi^-eo) 



cos^ I 
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Figure 1: Potential V_{x) (thick line) and the wave functions iPq{x), ^'^{x), 
'1^2 (solid line, short-dashed line and long-dashed line respectively) at the 
interval x G [0, 27r]. Here eo = 1, Cq = 0.05, Cf = 0.3, C2 = 1.3 are used. 

Potential V-{x) and the wave functions (x), ipi{^)y "^2 (^) W^- 
sented at the figure [TJ 

Because wave function ipQ^x) does not have nodes, eigenstate with en- 
ergy Eq = is a ground state of this potential. The wave functions ipi^x) 
and i^) have two nodes per interval of periodicity, then eigenstates with 
energies and E2 describe the limits of the second forbidden energy band. 

This quasi exactly solvable potential belongs to the class of QES poten- 
tials presented by Turbiner in his paper [S] in the following form 



in the case of?T, = l,a = l/2;aisa free parameter of quantum mechanics 
problem. 

Now let us consider supersymmetric partner of the potential V-{x): 



0} cos^(2aa;) — 2Q;a(2n + 1) cos(2aa;)). 



(78) 



^+(^) = ^0 + o^o - 1 - Vei^cosx - eoci cos^x + 



7 

^ ttj cos* X 

1=0 



(79) 
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2 ^ 6j cos* X 



i=0 
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ao = IQ44 

ai = -8y/e^eo{2 - 5eo + 2eg) 
a2 = -12eo(l-2eo-2e2 + 4ei]) 

as = 8v/i^(l + 3eo - I2el + 6^) , . 

a4 = l + 16eo-48e2(l-e2) ' 

as = -6v/i^(l + 2eo - I2el + Se;]) 
ae = -8e?eo(3 + 2eo) 

' bo = Seieg 
bi = 8el^/e^ 
62 = -2e2(l-12eo + 16e2) 
h = -8eov/eoer(3eo - 1) 
< 64 = eo(l + 10eo-48e2(l-eo)) • (81) 

65 = 2v/i^(l - 8eo + 12e2) 

66 = -2e?(-l-4eo + 16e2) 

67 = -8e?ye;5e7 
^ bs = 8eoef. 



Since we know eigenfunctions ipi (x) and ^/^^ (x) of Hamiltonian if_ , us- 
ing supersymmetric relations (|7j) we can find the wave functions ipt{x) and 
1(^2 (x), which are eigenfunctions of the Hamiltonian if+ with the correspond- 
ing energy values = eo and E2 = eo + Ci: 



COS' 




2 X] f^o^' ^ 



< 



i=0 

3 

^ nii cos* X 



(82) 



?/^^(x) 




— * 7—0 

2 smx^ 



2 ^ ni cos' a; 
1=0 



where 



ki = 4V2eoei /i = 2eo 



< A;2 = -y2(8eoei - 1) { k = 2{1 - Aeo)./!^ 



(83) 



fcg = -4V2eoei /s = 2ei 



^ k4 = 4A/2eoei [ ^4 = 4ei^eoei, 
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Figure 2: Potential l + (thick line) and the wave functions ipti^), i^ti^) 
(solid line and dashed line respectively) at the interval x G [0,27r]. Here 
eo = 1, = 0.2, = 0.7 are used. 



mo = -4V2eoei(4eo - l)(ei - y/e^) 

mi = 2v^(7?r - V^)(8eg - Meg + 7eo - 1) . . 

m2 = -v^(v^-ei)(l-4eo-4e2 + 16e3) ' ^ ^ 

ma = -4ely/2^{^ - y/^){4eo - 1) 

no = 2eoy/e^ 
rii = eo 

< na = -2(eo + ei)7i^^ (85) 

na = -ei 
^ n4 = 2ei^eoei. 

Thus we obtain QES potential VV(a;) ()79j) with two exactly know eigen- 
states El = cq, i'tix) and E2 = eo + £2, '?/'2'(^) given by (jH2I)- Potential 
V^(x) and the wave functions ipti^)^ i'ti^) presented at the figure |21 

Because the wave functions 'iptix) and V'^l^) have two nodes per period- 
icity interval, the eigenstates with energy values and E^ describe limits 
of the second forbidden energy band. Note, that QES potential (fT9|) does not 
belong to the general Turbiner's case and is completely new. 



20 



4 Conclusions 



In the present paper we have extended the SUSY method of constructing 
well-hke QES potentials with three known eigenstates potentials for the case 
of periodic potentials. 

Thus, periodic function U{x) generates quasi exactly solvable potential 
V-{x) with three known eigenstates iPq{x), ipi{x), '02 (^) quasi exactly 
solvable potential V+{x) with two known eigenstates 'ipi{x), ip2{x). 

Since we are interested in the real potential energy, condition R{x) > 
should be satisfied. To provide free of singularities potential energy and ex- 
tended bounded wave functions, generating function U{x) must have second 
order zero at the middle of the periodicity interval Xm and must be even func- 
tion with respect to this point. U{x) may have first-order zeros at the other 
points of the periodicity interval and should not have zeros of the highest 
order. The derivative of the U"{x) at the point Xm should satisfy the condi- 
tion U"{xm) = SeoCi. It is necessary to use solutions for the superpotentials, 
potentials and wave functions with opposite signs from the left and right 
sides with regard to point Xm to obtain continuous extended wave functions. 

As an example of the above described method starting from the gener- 
ating functions U{x) = 4eoei sin^ x we have obtained QES periodic potential 
V^{x) = eo — 1/2 + l/4(eoei — Gy^eoCi cosx — eoCi cos2x), which is trigono- 
metric extension of the well known Razavy QES potential, with three known 
eigenstates = 0, = eo and E2 = eo + ei, where ei = eo — 1/2 and 
eo is a free parameter. This potential belongs to the class of QES potentials 
presented by Turbiner at 0. Eigenstate with energy Eq = is the ground 
state of this potential. Eigenstates with energies E^ and E2 describes the 
limits of the second forbidden energy band. 

The supersymmetric partner V+{x) of potential V-{x) gives us a new 
QES periodic potential for which we know two eigenstates Ef = eo and 
E2 = €0 + 62 in the explicit form, where ei = eo — 1/2 and eo is a free 
parameter. This eigenstates describe the limits of the second forbidden energy 
band. 

Author is grateful to V. M. Tkachuk for enlightening suggestions, helpful 
comments and discussions. 
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